One- Way Quantum Computing in the Optical Frequency Comb 
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One-way quantum computing allows any quantum algorithm to be implemented easily using just 
measurements. The difficult part is creating the universal resource, a cluster state, on which the 
measurements are made. We propose a radically new approach: a scalable method that uses a 
single, multimode optical parametric oscillator (OPO). The method is very efficient and generates 
a continuous-variable cluster state, universal for quantum computation, with quantum information 
encoded in the quadratures of the optical frequency comb of the OPO. 
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Introduction Quantum computing (QC) is a fascinat- 
ing endeavor, ripe with promises of exponential speedup 
of particular mathematical problems such as quantum 
system simulation [1] and integer factoring [2]. Because 
implementing QC requires exquisite control of each single 
quantum memory unit (e.g., qubit) in a large-size regis- 
ter, practical QC is therefore faced with the daunting 
challenges of overcoming decoherence and achieving scal- 
ability [3] . Recently, the invention of one-way quantum 
computing introduced a new paradigm for quantum in- 
formation processing [4] based on teleportation alone [5] . 
In the traditional "circuit" QC model [6], physical quan- 
tum systems carry quantum information and undergo 
controlled unitary evolution; in the one-way QC model, 
quantum information exists virtually in a cluster state [7 
and is manipulated through a sequence of local measure- 
ments. The choice of measurement basis and the mea- 
surement results fully determine the quantum algorithm. 

The appeal of one-way QC is that it consolidates 
most of the challenging work into creating the univer- 
sal resource — the cluster state — and that it only re- 
quires local measurements. In addition, some one- 
way QC schemes admit very high fault-tolerance thresh- 
olds [8], and experimental realizations with four qubits 
have already been achieved [9, 10]. Efficient methods 
of creating large-scale cluster states are still needed, 
however, for practical implementation to be realistic. 
Here, we describe a radically new approach to scalabil- 
ity: a "top-down" method to produce large continuous- 
variable (CV) cluster states using a compact experimen- 
tal setup. The interest of CVs is their natural implemen- 
tation by squeezing (quantum noise reduction) in quan- 
tum optical systems. Photons are also less prone to de- 
coherence than, say, atoms, due to their lower propen- 
sity to interact with the environment. Several studies 
have established the use of photonic CVs for teleporta- 
tion [11, 12], QC [13], quantum error correction [14], clus- 
ter states [15, 16, 17, 18], and one-way QC [16]. 



CV cluster states from a single OPO Our implemen- 
tation of a scalable one-way quantum processor uses one 
multimode optical parametric oscillator (OPO). Initial 
proposals for constructing CV cluster states [15, 16] in- 
volved inline squeezers (seeded OPOs) [19], which are 
difficult to implement. A more viable method relics on 
the Bloch-Messiah decomposition [20] and uses N vac- 
uum squeezers followed by an C(7V 2 )-port interferome- 
ter [21]. Our method improves further, requiring only 
a single OPO and no interferometer [17], and provides 
huge scaling potential. 

The OPO combines two essential elements. The 
first one is an optical cavity (e.g., two facing mirrors) 
whose spectrum of resonant frequencies forms an op- 
tical frequency comb (OFC), so called because of the 
equal spacing between modes. Considered as a quan- 
tum system, the OFC is a large collection of indepen- 
dent quantum harmonic oscillator modes, or "qumodes" 
(a term used in analogy with "qubit"). Quantum infor- 
mation is encoded in the quadrature field variables of 
the OFC, which are analogs of position and momentum 
for a mechanical oscillator. The OFC can comprise mil- 
lions of modes and has outstanding classical coherence 
properties that have found groundbreaking applications 
in the revolutionary and now ubiquitous use of mode- 
and carrier-envelope-phase-locked femtosecond lasers in 
time/frequency metrology [22, 23]. 

The second crucial element of the OPO is the nondissi- 
pative nonlinear medium placed in the cavity. In a basic 
OPO, this medium is pumped by a monochromatic field 
and promotes downconvcrsion, i.e., the simultaneous ab- 
sorption of a pump photon at frequency lo p and emission 
of a photon pair at OFC frequencies u) m and uj n such 
that lu p = u> m + u> n , as well as upconversion, the reverse 
process. Such interactions yield bipartite CV entangle- 
ment of the OFC qumodes at frequencies u) m and uj n [24] . 
In the OPO we propose to use, the nonlinear medium is 
specifically engineered to quasiphasematch [25, 26] sev- 
eral such interactions simultaneously [27]. This, along 
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with a polychromatic pump, allows one to "write" an en- 
tangled network onto the OFC using pairwisc couplings. 
We show in this paper that this network can be made 
to precisely constitute a large CV cluster state, universal 
for one-way QC, and whose scaling requires no increase in 
pump complexity (number of frequencies) and only lin- 
ear increases in pump intensity and nonlinear gain band- 
width. 

A CV cluster state [15, 16] is defined as any mem- 
ber of a family of squeezed states indexed by an over- 
all squeezing parameter r > for which the variance of 
each component of (p — Aq) tends to as r — > oo [21]. 
Here, q = (t?i, . . . , <7at) t and p = (pi,... ,pn) T are vec- 
tors of amplitude and phase quadrature operators, re- 
spectively, and A is the (weighted, undirected) adjacency 
matrix of the cluster state's TV-node graph [17]. The 
infinite-squeezing limit is not achievable by any finite- 
energy state, an important point we address later. 

We consider an OPO that implements the following 
Hamiltonian, in the interaction picture and assuming a 
classical undepleted pump [17, 28]: 

H(A) = ihK22 A mn{aln a \i ~ Gro^n) , (1) 

where k > is an overall nonlinear coupling strength 
(squeezing parameter per unit time). The nodes of the 
graph described by A correspond to OPO qumodes and 
any (m,n) edge is weighted by A mn , whose magnitude is 
the qumode coupling strength in units of k and whose 
sign indicates downconversion if positive and upconver- 
sion if negative. In previous work [17], we showed that 
Ti.(A) always generates a CV cluster state whose graph is 
generally not given by A. If A is an orthogonal matrix, 
though, the graph is given by A — a fact we previously 
used to construct an H(A) that generates large sets of 
very small (2 x 2 or 2 x 3) CV cluster states [29]. Be- 
ing disconnected, these cluster states are not universal 
for QC, but requiring orthogonality is a useful simpli- 
fication. Here we construct an orthogonal A for which 
H(A) generates a QC-universal CV cluster state. 

Labeling rows and columns of A by sequential OFC 
modes, 7i(A) is easy to implement experimentally when 
A is a Hankel matrix (i.e., has constant skew diagonals). 
Any pump frequency u> p satisfying w p = uj m + u> n , to- 
gether with the assumption of a constant interaction 
strength, generates a constant skew diagonal in A and 
sets all its entries to the same value, fixed by the pump 
power. Additional pump frequencies generate additional 
skew diagonals, resulting in A having Hankel form. This 
connection, along with a useful shorthand for Hankel ma- 
trices, is illustrated in Fig. 1 and in the supplemental 
movie [31]. 

QC-universal CV cluster state We desire a CV clus- 
ter state that is universal for one-way QC [16], whose 
graph is bicolorable [17], and whose adjacency matrix is 
orthogonal and Hankel for experimental simplicity [29]. 
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FIG. 1: Hankel shorthand and pump specification. A Hankel 
matrix is uniquely specified by the entries along the top and 
down the right side. We collect these entries into a shorthand 
vector, representing the entire matrix itself, with the top-right 
entry set off with slashes. When A is Hankel, its shorthand 
vector immediately specifies the pump spectrum required to 
implement Ti.(A). Each nonzero entry in the shorthand vec- 
tor denotes the amplitude of a frequency in the pump, each of 
which generates CV entanglement between pairs of qumodes 
in the OFC symmetric about half that frequency. This ac- 
counts for all couplings prescribed by A. 



Orthogonality of the adjacency matrix (AA T = 1) for an 
undirected graph (A = A T ) yields A 2 = 1, or 

'A ;! K, X-F/.W 6 jk . (2) 

i 

Eq. (2) has a geometric interpretation: (A™)^ represents 
the sum of the weights of all n-length paths from node j 
to node k, where the weight of such an "n-path" equals 
the product of all edge weights along the path. Eq. (2) 
requires that all 2-paths that begin and end on the same 
node have weights that sum to 1, while those that link 
different nodes have weights that cancel out. 

The usual candidate for a QC-universal graph is the 
square lattice [4]. This graph is irregular at the bound- 
aries, so Eq. (2) cannot be satisfied for any real-valued 
weighting of it. Applying toroidal boundary conditions 
will make it regular, but orthogonality is still prohibited 
because there exist pairs of distinct nodes connected by 
only one 2-path, for which the sum in Eq. (2) collapses 
to a single term (which must be nonzero). 

While real-valued weights cannot satisfy Eq. (2) for 
a toroidal square lattice, matrix-valued weights can. In 
such a case, Eq. (2) becomes J^i ^ji^ik = <5jfcli where the 
"entries" A,,-; are themselves mxm matrices. This means 
A is now an adjacency matrix on an {mN)-node graph. 
On the other hand, treating the mxm blocks as single 
entries, A is the matrix-weighted adjacency matrix for 
what we call a supergraph, which has macronodes con- 
sisting of m individual nodes each. Fig. 2 illustrates the 
meaning of matrix-weighted edges between macronodes 
for the simple case of a ring supergraph, where m = 2. 

We henceforth promote the toroidal square lattice to a 
supergraph, illustrated in Fig. 3. The lattice has degree 
four, so we choose our matrix-valued weights to be four 
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FIG. 2: Matrix-valued weights and supergraphs. The matrix- 
valued weights 7r + = 5 ( + + ) and tt~ — | ( 1 7_ ) , where ± 
stands for ±1, connect the macronodes of a "ring" supergraph 
(left). The entries in n ± specify the real- valued weights in the 
actual "crown" graph (right) that connects the underlying 
physical nodes. Measuring q for each of the physical nodes 
in the top layer of the crown leaves the bottom layer in a 
uniformly weighted ring-graph CV cluster state. 
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where ± stands for ±1. Note that the toroidal embedding 
in Fig. 3 involves a twist of one unit in both dimensions, 
a subtlety of little physical consequence but which allows 
for a Hankel A at the supergraph level (using the Hankcl 
shorthand defined in Fig. 1): 

a = [o^n^o^n^o^n^o/n^o^n^o^n^o^-n^o] , (4) 

where k is a list of k zero-matrices, the size of which 
is governed by context (4 x 4 in this case), there are 
M 2 macronodes (with M even), u = M — 1, and v = 
M 2 - 2AI - 3. (The negative sign on II 3 is used, with 
foresight, only to make the node renumbering below work 
out properly.) Orthogonality (AA T = 1) is easily verified. 

The M-indexed family of these cluster states is uni- 
versal for CV one-way QC. To see this, first measure q 
on three physical nodes per macronode to reduce the su- 
pergraph to a uniformly-weighted graph with the same 
structure (see Fig. 3) . Then cut open the toroidal lattice 
by measuring q along two orthogonal meridians to form 
an ordinary lattice, known to be universal [16]. 

The matrix A, while Hankel at the supergraph level, is 
block-Hankel (with 4x4 blocks) at physical node level. 
A suitable renumbering of nodes will reduce the blocks 
to 2 x 2 while maintaining the block-Hankel structure: 

A = [O s ,7r~,O*7r + ,O s ,7r + ,0,7r~,O s ,7r~ O*,7r + ,O s ,7r",0/7r + / 

0>- 0*7r + ,0> + ,0,7r- OV" 0*7r+OV",0] , (5) 

where 7r ± are defined in Fig. 2, s = 2AI — 1, t = 
M 2 — AM — 3, and = indicates equality up to node renum- 
bering. These cluster states are QC-universal, bicol- 
orable, and orthogonal — but still only block-Hankel. This 
is nevertheless sufficient for simple experimental imple- 
mentation, using a method we now describe. 




FIG. 3: Toroidal lattice supergraph and underlying graph 
structure. Each macronode in the supergraph (left) consists 
of four physical nodes, and each color corresponds to one of 
the four matrix- valued weights IP from Eq. (3). Entries in 
IF specify the real-valued weights connecting the underlying 
physical nodes (right). Measuring q on each physical node 
in three of the four "layers" leaves the remaining layer in a 
uniformly weighted QC-universal toroidal lattice cluster state. 



Experimental implementation In Eq. (5), each ■n ± 
block corresponds to a single pump frequency. Such 
couplings can be implemented using the two orthog- 
onal polarizations of an optical field at a given fre- 
quency, as was experimentally demonstrated by simulta- 
neously quasiphasematching polarization-sensitive inter- 
actions ZZZ, ZYY, and YZY/YYZ (first letter is pump 
polarization) in a periodically poled KTiOP04 (KTP) 
crystal [27]. The difference between ir~ and 7r + is thus 
a 180° phase-shift in the Y-polarized pump mode. A 
narrowband pump polarized at ±45° in the (ZY) plane 
implements a it^ skew-diagonal band in A. Equation (5) 
therefore translates into a single OPO pumped by exactly 
15 frequencies. Interactions with OFC modes outside the 
desired subset must also be strictly suppressed, which can 
be done by cavity mirror design and/or by quasiphase- 
matching bandwidth design [20]. While the pump spec- 
trum is relatively complex — requiring 15 frequencies — 
that number is constant with respect to the lattice size, 
making this construction extremely scalable. 

Note that, in principle, additional physical parameters 
like wave vector direction or transverse mode structure 
could be used to directly implement block-Hankel A's 
with larger blocks (e.g., 4x4). This would reduce pump 
complexity but require a more sophisticated OPO. 

Finally, we provide realistic estimates for the scaling 
potential of the CV cluster state with N macronodes 
(N = M 2 ) and constant overall coupling strength k. The 
pump spectrum complexity and the OFC are indepen- 
dent of N. Only two quantities scale linearly with N: the 
overall pump power and the bandwidth of the nonlinear 
coupling. The number of pump photons must increase 
with the number of entangled qumode pairs (number of 
graph edges), which grows linearly with N in a square 
lattice. Two-mode CV entanglement can be obtained at 
a few milliwatts pump power, with the upper limit be- 
ing the optical damage limit, which in KTP is at least 
several watts (focused, continuous- wave), thus yielding 3 
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orders of magnitude of scaling range. The bandwidth of 
the nonlinear coupling (100 GHz to 1 THz) must encom- 
pass the whole desired set of qumodes, separated by the 
cavity's free spectral range (100 MHz for a 1.5-meter- 
long cavity), which yields 3 or 4 orders of magnitude. 
These figures reflect ordinary — rather than state-of-the- 
art — performance and do not account for other interest- 
ing avenues such as implementing nonlinear couplings in 
slow light media. These estimates indicate our approach 
has a quite realistic potential for significant scaling. 

Finite squeezing and CV fault tolerance The finite- 
squeezing approximation is a special case of more gen- 
eral considerations of error correction and fault tolerance 
for one-way QC using CV cluster states. Certainly, more 
squeezing is preferable to less, but the amount required 
for any particular QC task remains an open question. 
Consequently, it's unclear how the nonlinear coupling 
strength k in Eq. (1) will need to scale with N for any 
particular QC application. Our results are nonetheless 
compelling because as N increases, the new interactions 
generated have the same squeezing strength as existing 
ones — i.e., existing squeezing is not "redistributed" to 
the new pairings as the lattice and pump power grow. 
Moreover, scalability does enable quantum encoding re- 
dundancy. While finite squeezing errors can be mitigated 
in medium-sized proof-of-principlc experiments with CV 
cluster states [16], and some work also addresses CV er- 
ror correction in general [14, 30], the task of rigorously 
establishing a fault tolerance threshold for CV one-way 
QC is an important open problem. We hope to spur fur- 
ther investigations along these lines. 

Conclusion We have presented a theoretical break- 
through that opens the door to large-scale experimental 
generation of a universal one-way quantum computing re- 
source, using existing technology The entangled states 
produced by this method will also be objects of interest 
in the study of entanglement at mesoscopic scales. Fur- 
ther simplification of the method may be possible using 
additional degrees of freedom, such as spatial modes. 

This method of entangling — in one fell swoop — a large 
optical frequency comb into a continuous-variable clus- 
ter state is the first "top-down" approach proposed for 
one-way quantum computation using optical encodings 
of information in either discrete or continuous variables. 
While open questions remain about the effects of finite 
squeezing on scalability for particular quantum comput- 
ing tasks, the unprecedented scalability of this method 
encourages further theoretical research. Experimental 
implementation is already underway. 
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